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First results with non-perturbative fermion improvement* 
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We present initial results for light hadron masses and nucleon structure functions using a recent proposal for 
eliminating all 0(a) effects from Wilson fermion simulations in the quenched approximation. With initially limited 
statistics, we find a much more linear APE plot and a value of the axial coupling qa nearer to the experimental 
point than with comparable runs using unimproved Wilson fermions. 



1. NON-PERTURBATIVE 
IMPROVEMENT 

It is now common to seek to improve Wilson 
fermions by the addition of a term of the form 



Ssw = ^/tffcswa/ ^(x)a flv F flv (x)ip(x] 



(1) 



in order to reduce O(a) lattice cut-off effects 
present in the action. The original proposal jj] 
was to use csw = L which corresponds to the 
tree-level perturbation theory result; this elmi- 
nates terms up to 0(ag 2 ). The elimination of 
tadpole contributions requires a value of csw 
some 50% larger depending on the coupling |Q, 
though the leading corrections here remain for- 
mally 0(ag 2 ). 

It has recently been suggested || that csw can 
be chosen via some suitable physical condition so 
as to remove all 0(a) effects. Here we present ini- 
tial results using this action as part of our project 
for a non-perturbative calculation of nucleon ma- 
trix elements Hi. 



2. IMPLEMENTATION 

Our initial runs used a 16 3 x 32 lattice at [3 = 
6.0 with csw = 1.769 ||. Our implementation 
runs on a Quadrics QH2 parallel computer with 
an 8 x 8 x 4 topology. The improvement term 
in equation (Q) appears in the site-diagonal part 
of the action; the major overhead in our case is 
multiplication by this term during inversion of the 
fermion mass matrix. In our gamma matrix basis, 
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*Talk presented by P.W. Stephenson at Lat96, St. Louis. 



(where A, B are 6x6 matrices, i.e. two-spinors 
with colour) so that instead of a 12 x 12 multipli- 
cation we have two 6x6 multiplications and two 
inexpensive co-ordinate transformations. This re- 
duces the overhead for the improvement in the 
inverter from 45% to 30%. Also, the inverse of 
the term in equation (3) is required on half the 
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lattice due to the red-black preconditioning; we 
now have to invert two 6x6 instead of a 12 x 12 
matrix. However, this is only required once for 
each propagator inversion, so is a less significant 
factor. 

At (3 = 5.7 we see large numbers of configura- 
tions where the fermion matrix inversion fails to 
converge even after 2000 minimal residual sweeps. 
For csw > 1 these amount to more than 5% of 
the total. For j3 < 6.0, no non-perturbatively 
calculated csw is available H for essentially this 
reason. 

3. HADRON MASSES 




1.3 - 



■ 

1.2 - 



' 0.0 0.2 0.4 0.6 0.8 1.0 

( m „ /m p) 2 

Fi gure 1 . The APE plot for the unimproved (unfilled 
symbols) and cgw = 1.769 improved (filled) data. The 
lattices used are 16 3 X 32 and 24 3 X 32 at j3 = 6.0. The 
physical and heavy quark points are shown as squares. 

We have calculated pion, rho and nucleon 
masses for k = 0.1324, 0.1333 and 0.1342 from 
125 configurations, and on a 24 3 x 32 lattice at 
0.1342, 0.1346 and 0.1348 with similar statistics. 
The APE plot is shown in figure |], where we com- 
pare the results with previous runs using csw = 
at k — 0.1515, 0.1530 and 0.1550 with large statis- 
tics (some 5000 data points) on the smaller lattice 



and 0.1550, 0.1558 and 0.1563 again with around 
100 configurations on the larger lattice (the high 
statistics and light quark results are both previ- 
ously unpublished) . The k's have been chosen to 
lead to similar pion masses in the improved and 
unimproved data sets; we also find that the nu- 
cleon masses are similar. However, the rho mass 
with improved fermions is some 20% larger for 
the higher mass points, which are thus shifted to- 
wards the origin in comparison with the Wilson 
case. 

As may be seen, we find that the unimproved 
Wilson hadron masses eventually converge with 
the improved values, and indeed the extrapo- 
lated values are comparable, although the extrap- 
olation of the ratios deviates considerably more 
from the linear in the former case. We have also 
found extrapolation of the individual masses to 
the chiral limit to be more linear for the improved 
fermions. 

4. STRUCTURE FUNCTIONS 

For a full non-perturbative calculation of the 
structure functions we need two additional ingre- 
dients. 

First, we require operator improvement which 
needs to be calculated separately for each opera- 
tor. The perturbative case has 

V = iyyui> — ► ipluip 

(4) 

+ bvai-tpD^tp - <3„V<vV0> 

A = ipj^ 5 4> -> V'7m75?/' 

_ - (5) 
+ b A a(dfj,ipj 5 ip + ipcrp V D„75i>); 

in this preliminary calculation we have set by = 
b A = 1/2, corresponding to the 'rotation' valid in 
the csw = 1 case. Non-perturbatively, we need to 
determine the coefficients of each term separately. 

Secondly, we also need to determine the renor- 
malisation constants Zo non-perturbatively. At 
present we use a perturbative calculation assum- 
ing the by = bA = 1/2 rotation but with csw = 
1.769. A point to note is that definitions of Zo 
vary, depending on whether the residual effect of 
the rotation in the chiral limit is absorbed into Zo 
or not: we follow Martinelli et al. ||. This is dif- 
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ferent to the prescription being used by the AL- 
PHA collaboration for these calculations, hence 
their results presented at this conference |(| are 
not immediately applicable to us. The calcula- 
tion of the Zo is then similar to that presented 
in reference 0. 

A fully non-perturbative calculation of the 
renormalisation constants is currently in progress. 

4.1. Local vector current 

We have calculated the local vector current and 
extrapolated to the chiral limit, where the per- 
turbative value for Zy was calculated. As this 
current is conserved in the continuum we require 
only to recover the values 2 for the up quark and 
1 for the down quark content of the proton. This 
acts as a test of the consistency of our procedure. 
The results are 2.0(3) and 1.0(2), which are en- 
tirely satisfactory. 

4.2. Axial vector current 

The axial vector current Aq is shown in fig- 
ure 0. This was formerly (in the naive quark 
model) connected with the spin contribution of 
the quarks, though emphasis is now on the viola- 
tion of the OZI rule implied by the low values for 
the flavour singlet operator M. 

Our main interest is in the combination gA = 
Au — Ad (the axial coupling) where disconnected 
fermion loops which we have not calculated cancel 
out. The experimental value is 1.26. With the im- 
provement, our calculated quantity has changed 
from 1.07(9) to 1.22(14). This is better but suf- 
fers from large errors. 

5. CONCLUSIONS 

We have performed a first QCD calculation 
with Wilson fermions improved up to 0(a 2 ). 
Hadron masses with improved fermions at (3 = 
6.0 extrapolate more linearly to the chiral limit 
than with ordinary Wilson fermions. 

Our preliminary analysis of the local vector 
current in the nucleon shows the procedure to 
be consistent, and our results for the axial vec- 
tor indicate a promising trend in our value for 
the axial coupling despite low statistics. A fully 
non-perturbative analysis remains to be done. 

The work was done on the APE Quadrics com- 
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Figure 2. The axial vector currents for up and down 
quarks, Am and Aq!, for both improved (A) and unim- 
proved (o) fermions, including the factor calculated in 
the chiral limit (solid points). 

puters at DESY-IfH, Zeuthen, Germany. We 
thank the computer center for technical help. 
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